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We investigate the capability of LISA to measure the sky position of equal-mass, nonspinning 
black hole binaries, including for the first: time the entire inspiral-merger-ringdown signal, the effect 
of the LISA orbits, and the complete three-channel LISA response. For an ensemble of systems 
near the peak of USA’s sensitivity band, with total rest mass of 2 x 10 6 M© at a redshift of ,7 = 1 
with random orientations arid sky positions, we find median sky localization errors of approximately 
~ 3 arcminutes. This is comparable to the field of view of powerful electromagnetic telescopes, 
such as the James Webb Space Telescope, that could be used to search for electromagnetic signals 
associated with merging black holes. We investigate the way in which parameter errors decrease 
with measurement time, focusing specifically on the additional information provided during the 
merger-ringdown segment of the signal. We find that this information improves all parameter 
estimates directly, rather than through diminishing correlations with any subset of well-determined 
parameters. 
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I. INTRODUCTION 

Among the gravitational-wave (GW) sources which the 
proposed Laser Interferometer Space Antenna (LISA) 
would observe, massive black hole binaries (MBHBs) 
stand out in several ways. They have the highest typical 
signal-to-noise ratios (SNRs) and will be observable out 
to redshifts of z ~ 10 or greater [1 3], making MBHBs 
the most distant sources that LISA will observe. Obser- 
vations of MBHBs figure prominently in LISA’s science 
goals [4] and may provide particularly rich opportunities 
for multi-messenger observation in conjunction with elec- 
tromagnetic instruments. The potential for such obser- 
vations depends particularly on LISA's ability to localize 
the gravitational-wave sources on the sky and thereby 
facilitate simultaneous or follow-up electromagnetic ob- 
servation [5]. 

LISA is an all-sky instrument in the sense that, at any 
one instant, it is sensitive to sources at most points on 
the sky. While this is advantageous from the perspec- 
tive of detection of unknown sources, it is not helpful 
for locating these sources on the sky. To build localiza- 
tion information, LISA takes advantage of the fact that 
MBHBs are observable for months or even years prior to 
merger. During this time, LISA’s orbit around the Sun 
introduces both frequency and amplitude modulations to 
the signal that can be used to determine source location 
[6]. The majority of studies of LISA's source location 
abilities [7—13] focus on this mechanism and neglect the 
final merger-ringdown portion of the MBHB waveform. 
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which is difficult to model and too brief to experience any 
modulation from LISA’s orbit. 

Recent work [14, 15] has shown that the short-lived, 
but stronger and higher-frequency, merger-ringdown por- 
tion of the MBHB signal alone can achieve compara- 
ble sky localization to the long inspiral portion without 
any additional information from LISA’s orbital modula- 
tion. The mechanism for extracting this information is 
believed to be related to the complex response of LISA 
to gravitational waves at the upper end of the measure- 
ment band (above the inverse light crossing time or about 
60 mHz). In this paper, we investigate the combined 
effect of both localization mechanisms. Specifically, we 
expand our previous analysis to include the full inspiral 
signal, including orbital effects. We also include a third 
signal channel, T, relevant only at higher frequencies, 
which was excluded from the earlier analysis for techni- 
cal reasons. 

The outline of the paper is as follows: In Section II, 
we review our methodology, which closely follows [15]. 
Section II A describes the construction of the inspiral- 
merger-ringdown (IMR) waveform using the “IRS-EOB" 
model. Section II B describes our model for the LISA 
instrument, including the response function and noise. 
Finally, Section II C explains how we calculate parameter 
errors using the standard Fisher matrix formalism. 

In section III, we examine LISA’s ability to localize 
sources, including the improvements that are gained as 
different levels of realism are added. We find that for an 
ensemble of systems with total rest mass of 2 x IQ 6 M© at 
a redshift of z = 1 , observing the merger-ringdown (MR) 
signal with a stationary detector provides better localiza- 
tion information than observing the inspiral signal with 
an orbiting detector. We then consider the observation 
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of a complete IMR signal with an orbiting detector and 
find the final error to be only a negligible improvement 
on the MR signal, stationary detector result. Measur- 
ing the long-duration inspiral is essentially irrelevant to 
achieving LISA’s localization potential. Nevertheless, the 
orbital information contained in the inspiral is still im- 
portant for localizing signals in advance of merger, which 
may be critical for performing counterpart searches. We 
find, in agreement with previous results, that LISA can 
localize most MBHBs to within the Large Synoptic Sur- 
vey Telescope (LSST) field of view (FOV) several months 
before merger. With the complete IMR signal, ~ 50% 
of MBHBs in our ensemble can be localized into the 
much smaller FOV of the James Webb Space Telescope 
(JWST). This achievement is possible mainly because of 
the extra information provided by the T channel. 

In Section IV, we demonstrate that the rnerger- 
ringdown portion of the waveform provides additional 
information beyond the inspiral for all parameters (al- 
though the improvement in mass measurements is rela- 
tively small). Furthermore, we show explicitly that the 
reduction in sky localization errors is not primarily due 
to improved measurement of coalescence time. We sum- 
marize and conclude in Section V, focusing on how our 
results may impact alternative designs for LISA. 


II. METHODOLOGY 


Working in geometrized units where G = c — M = 
mi + m -2 = 1, the source model depends only on the 
mass ratio q = rn [ / m -2 (where m \ < roo) and the spins. 
We generate a waveform time series using the IRS-EOB 
model with a signal cadence corresponding to a quarter 
wavelength at the highest frequency reached by the £ — 4, 
m = +4 harmonics. After the source calculation, we 
transform the waveform from the source coordinates to 
the solar system barycenter (SSB). The SSB waveform 
h b is given by: 
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where ~ 2 Y( m are the spherical harmonics of spin-weight 
—2 [18] and ht m is the dimensionless strain decomposed 
in that basis in the source frame. This transformation 
depends on six parameters: the redshifted total system 
mass M = M 0 ( 1 + z), the luminosity distance £>/,, the 
coalescence time t c , and three angles describing the orien- 
tation of the binary. The coalescence time, t c , is defined 
as the time at which the peak of the waveform reaches 
the solar system barycenter. For the angles, we use the 
inclination i (with the convention that /, = 0 corresponds 
to the line of sight being coincident with the orbital axis 
of the binary), initial orbital phase <j> 0 , and polarization 
phase 


The procedure used to estimate the errors with which 
LISA will measure astrophysical source parameters of 
MBHBs closely follows that of our previous studies 
[14, 15]. We will provide a brief overview; those inter- 
ested in more detail should consult [15]. To assess pa- 
rameter accuracy, we require a model for the emitted 
waveform, a model for the detector response and noise, 
and a method for converting these pieces of information 
into a theoretical limit on the achievable accuracy in mea- 
suring parameters with LISA. 


A. Waveform Model 

We use a complete IMR waveform model [16] tuned to 
match the available numerical simulations for nonspin- 
ning black hole binaries. This model, referred to as the 
IRS-EOB model, uses a conventional effective-one-bodv 
(EOB) Hamiltonian formalism to model the adiabatic in- 
spiral [17]. The merger-ringdown is modeled using the 
“implicit rotating source” (IRS) formalism, which treats 
the radiation source as a shrinking rigid rotator. 

The waveform amplitude is calculated using a flux 
model that is constrained both to be consistent with the 
inspiral flux through 3.5 post-Newtonian (PN) order and 
also to vanish as it approaches the ringdown frequency 
(referred to as “Model 2” and given by Eq. 19 in [16]). 
The waveform phase is fit to a physically motivated func- 
tional form (see Eq. 9 in [16]). 


We focus on a specific case that has appeared elsewhere 
in the literature (e. g. [7, 9]): an equal-mass, nonspinning 
MBIIB with rest mass M a = 2 x 10 6 M® at 2 = 1, 
observed for the final year prior to merger. This mass is 
interesting because it lies near the most sensitive range 
of the LISA instrument [3], merger-tree models predict 
LISA will observe a significant number of events in this 
range [19], and excellent observational evidence exists for 
black holes of this size [20 22] . 

The sky location of the binary, described by the ecliptic 
latitude j3 and longitude A, is applied by the instrument 
response, which we discuss below. We use the vector 
A“ = (In M, In £)/_,, 3, A, t, d> 0 , ip, t c ) to denote the com- 
plete set of variable parameters. We exclude the mass 
ratio q from A'% as it is not varied when computing pa- 
rameter uncertainties, a procedure consistent with that 
used in [14, 15, 23]. 


B. Instrument Model 

The instrument model includes the transformation of 
the SSB waveform Iib to the LISA instrument outputs as 
well as a model of the instrument noise. The instrument 
outputs are the {A, E, T] TDI combinations described in 
[15] and based on the orthogonal {A, E, T} combination 
[24! . The response is computed by employing Synthetic 
LISA [25] to generate the {X, Y. Z} TDI combinations 
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[26] and then converting to {A, E, T} using 
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We use an analytic model of the noise in the {*4, E, T} 
TDI channels, obtained following the procedure in [27]. 
The expressions for the one-sided spectral densities are 


(2) 


s a,e = 2 sin 2 (4>) [2 (3 + 2 cos(<3?) + eos(2 #)) S pm + (2 + cos(4>)) 5 op ] , 

S f = 8sin 2 ($)sin 2 («h/2) [4sin 2 ($/2)5 pm + 5 op ] , (3) 
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where <I> = 2-irfL/c, L is the arm length, and S pm and 
,S’ 0 p are the one-sided spectral densities of the proof-mass 
acceleration and optical path-length noises, 


Spm = 2.5 x 10- 
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In our prior studies [14, 15], we neglected the T chan- 
nel due to inconsistencies that were observed between the 
analytic noise model for Sf in (3) and the noise spectrum 
obtained from a time-domain simulation of Sf performed 
with Synthetic LISA. We have now verified that the 
two models are consistent for / > 0.5 mHz. We now 
include contributions from the T channel at these fre- 
quencies where noted. We expect this to be a reasonably 
accurate approach since the response of the T channel 
vanishes at low frequencies. 

For the foreground of gravitational waves from unre- 
solved compact binaries, we use the model developed in 
[28], by adding 

Sga] = [4$ sin($)] 2 S con f (5) 

to the expressions for S& g. Sconf is a signal confusion 
noise estimate for a related measurement variable, taken 
from Eq. 14 of [28]. No gravitational- wave foreground is 
added to the T channel since S con { ■C Sf at the frequen- 
cies where T is included. 


C. Error estimation 

To approximate the achievable measurement accuracy 
for LISA, we apply the Fisher matrix formalism. While 
not as robust or informative as more sophisticated meth- 
ods such as Markov Chain Monte Carlo [29] , the compu- 
tational efficiency of the Fisher matrix approach is well 
suited for this type of performance estimation work. Po- 
tential pitfalls have been discussed in the literature [30]. 


The Fisher information matrix is defined as 


r a & 
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where h is a channel of TDI strain, a and b are parameter 
indices, and (...[...) denotes a frequency-domain inner 
product weighted by the inverse noise spectral density 
S n ( f) of the TDI output n = A, E, T: 
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Here /i ow is the low-frequency cutoff, chosen specifically 
for each estimate (see Section III), and /high is the max- 
imum frequency present in our waveforms, as described 
in Section II A. The individual Fisher matrices for each 
TDI channel to be included in the analysis are summed to 
form the total Fisher matrix. Throughout this work, we 
calculate the parameter derivatives in (6) using second- 
order centered differencing (rather than one-sided as in 
[15]), with a fractional step size e a = AA a . We set 
A = l(r 4 for the coalescence time and A = 1()~ 6 for 
all other parameters. 

The covariance matrix, E ah , is well approximated by 
the inverse of the Fisher matrix for signals w r ith large 
SNR: 

E ab = (F° 6 ) ~ 1 [l + C7(SNR _1 )] , (8) 


so that a a = y S tt “ is the standard deviation of parameter 
A“. Because we use In M and In Di as parameters, the 
resulting uncertainties are fractional: 

er ln M = cr M /M , 

a \nD L = (T Dl/ Dl (9) 

The uncertainty of the other dimensionful parameter, t c , 
is not expressed fractionally, as one is generally interested 
in the absolute timing error, but we emphasize that one 
is free to make either choice. We note that a poor choice 
of units for t c could cause the Fisher matrix to have a 
large dynamic range and be ill-conditioned for inversion. 
We find that expressing t c in seconds generally yields 
computationally invertible Fisher matrices. 



4 


III. SKY LOCALIZATION 

As mentioned in Section I, the localization of a par- 
ticular MBHB source on the sky is critical for enabling 
coordinated observations with electromagnetic observa- 
tories. One goal might be to locate the host galaxy in 
which the MBHB merger occurred. If the galaxy can 
be identified and the redshift obtained, it can be com- 
bined with the luminosity distance determined from the 
gravitational waveform for use as a cosmological probe 
[5, 31]. The simplest way to identify a host galaxy would 
be for LISA to generate an error volume that only con- 
tained a single galaxy. Unfortunately, even the most op- 
timistic estimates of LISA’s localization ability suggest 
that this will not generally be possible [5]. Identifying 
the host galaxy within LISA’s error volume will require 
other information, perhaps the detection of a transient 
electromagnetic signal associated with the merger itself. 
The nature of this possible signal is a subject of active 
research [32, 33], and it is unclear if it exists, when it 
might peak relative to the GW signal, and what its spec- 
tral content will be. 

The observational effort needed to search for an EM 
counterpart to a MBHB merger will depend on the ra- 
tio between LISA’s sky localization error box and the 
field of view (FOV) of the EM asset used to perform the 
search. Roughly speaking, this ratio will determine the 
number of telescope pointings that are needed to cover 
the LISA error box. While some counterpart scenarios 
involve afterglows that persist long after coalescence, oth- 
ers involve signals that peak at or before coalescence and 
may fade quickly. For these signals, it is important to 
have some advance warning as to where on the sky (and 
when) the merger will occur, allowing for appropriate 
EM resources to be scheduled and pointed in advance. 
The sky localization information produced by LISA’s or- 
bit accumulates with time over the months or years that 
the MBHB is observed, allowing for such advance notices 
to be issued. 

Figure 1 summarizes various estimates for the devel- 
opment of precision in ecliptic latitude over time for a 
LISA observation of two equal-mass, nonspinning black 
holes with a combined rest mass M 0 = 2 x 10 6 M@ at 
2 = 1. All of the curves show the median results from 
an ensemble of 100 systems with random sky locations 
and orientations. Except where noted, the low-frequency 
limit of the inner product (7), f\ ow , was set to 0.03 mHz. 
The precise placement of this cutoff and the model of the 
instrument noise between the cutoff and the lower limit 
of the “official” LISA measurement band at 0.1 mHz can 
have a significant effect on the results [5, 9]. The curves 
are all shown as functions of time before t c , defined as 
the peak of the waveform. In addition, we also show to 
the right the final error achieved with the complete signal 
including ringdown. 

The purple dot-dashed curve utilizes PN waveforms 
and a simplified model of the LISA instrument response 
known as the low-frequency approximation (LEA). This 


approximation produces two channels of information 
which can be compared to the A and E TDI channels dis- 
cussed in Section II B. The estimate is generated by the 
same code used in [11, 12], with the exception that here 
we ignore the effects of spin precession. We include it pri- 
marily for comparison with our own method described iii 
Section II. Due to the breakdown of the post-Newtonian 
approximation, the evolution of' sky localization infor- 
mation is truncated at the time corresponding to the 
Schwarzschild innermost stable circular orbit (ISGO) or 
r = 6M. Therefore, this estimate is sensitive to sky local- 
ization information encoded in the waveform by LISA’s 
orbital motion but not sensitive to localization informa- 
tion produced by the interaction of the MR waveform 
with the high-frequency response of the instrument. 

The red dotted curve, based on our previous work [15], 
uses the IRS-EOB waveform model, as described in Sec- 
tion II A, and the Synthetic LISA instrument response 
described in Section II B. Due to the aforementioned 
problem (now solved) of spurious information in the T 
channel, only the A and £ TDI channels are included. In 
addition, the instrument response uses a stationary, non- 
orbiting constellation lying in the ecliptic plane instead 
of the usual LISA orbits. This estimate is therefore only 
sensitive to the sky localization information provided by 
the interaction of the MR signal with the high-frequency 
response of the detector. Lacking the orbital contribu- 
tion, we can safely remove a large portion of the inspiral 
by setting the low-frequency limit to f\ ow = 0.1 mHz for 
this single estimate only. This helps reduce computation 
time. 

For times before the ISCO, the errors in the orbits-only 
estimate (purple, dot-dashed) are far less than those in 
the stationary-LISA estimate (red, dotted), suggesting 
that the orbital mechanism for sky localization is the 
dominant one over these time periods. The dotted curve 
only “catches up” to the dot-dashed curve at ~ 100 sec- 
onds before the merger. However, by the time the full 
signal is considered, the stationary-LISA errors have be- 
come ~ 4 times smaller than the orbits-only errors at 
ISCO. In other words, for nonspinning, equal-mass bina- 
ries in this mass range, observing the MR signal alone 
gives better localization information than observing the 
entire inspiral. 

The remaining curves are new results, estimates which 
include both the complete IMR (IRS-EOB) waveform 
and a realistic, orbiting (Synthetic LISA) model of 
the instrument. In addition, the low frequency cutoff has 
been restored to 0.03 mHz in order to take advantage 
of the orbital modulation. The green dashed curve in- 
cludes only the A and E TDI channels, while the blue 
solid curve includes the T channel as well. 

The A , E estimate (green, dashed) tracks the orbits- 
only estimate (purple, dot-dashed) at times well before 
merger, thus confirming that our model is consistent with 
that of Lang and Hughes [11, 12] in the inspiral regime. 
It also tracks the stationary-LISA estimate (red, dotted) 
as the MBHB approaches merger. When the entire wave- 
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FIG. 1: Comparison of latitude uncertainty estimates for an equal-mass, nonspinning MBHB with total mass 2 x 1() 6 M@ 
and redshift z — 1. The purple dot-dashed estimate (based on [11, 12]) uses a post-Newtonian inspiral (I) waveform and a 
low-frequency approximation (LFA) to the detector response that includes USA’s orbit. The red dotted curve (based on [15]) 
uses the complete inspiral-merger-ringdown (IMR) waveform and Synthetic; LISA (SL) for the response function but neglects 
LISA’s orbit, all frequencies below 0.1 mHz, and the T channel. The green dashed curve uses the IMR waveform, the SL 
response, and orbital effects, but neglects the T channel. The blue solid curve includes all effects; the IMR waveform, the SL 
response with orbits, and all three TDI channels. The FOVs of LSST and JWST, as well as the location of the Schwarzschild 
ISCO, are shown for reference. 


form is used, these two estimates differ by only ~ 10%. 
The dashed curve contains both mechanisms for localiz- 
ing a source (orbits and high-frequency detector response 
to merger-ringdown), but it ultimately does not do signif- 
icantly better than the stationary-LISA response alone. 
Our ability to localize sources well in advance of merger 
depends critically on the orbital information. However, 
for the final achievable sky location accuracy in these sys- 
tems, the inspiral is essentially irrelevant and only the 
AIR signal need be taken into account. This result has a 
significant impact on potential new LISA designs, which 
may feature different orbits and shorter arms (and thus 
a higher frequency cutoff). We discuss these impacts in 
more detail in section V. 

The estimate with all three TDI channels (blue, solid) 
begins to diverge from the orbits-only estimate (purple, 
dot-dashed) at earlier times than the A , E estimate and 
reaches a value at coalescence that is about half that of 
the A , E estimate. This represents the additional con- 
tribution of the T channel, which is concentrated at high 
frequencies. It is neglected implicitly in the orbits-only 
estimate through the use of the LFA and explicitly in 
the other two estimates. While all three TDI channels 
are available in a fully operational, 6-link LISA, only two 
channels are available if one link is down. Using the green 
dashed line as a proxy for a 5-link LISA and comparing 
to the blue solid line, we can see roughly how much the 
sky position determination is affected by a non-optimal 
detector. 

Figure 1 also includes the FOV of two representative 
EM instruments that could be used for locating coun- 


terparts. The Large Synoptic Survey Telescope (LSST) 
[35] has a FOV of 9.6 deg*, which we have modeled as a 
square aperture ~ 186 arcmin on a side, and is designed 
to rapidly and repeatedly scan large areas of the sky. 
The figure demonstrates that LISA will localize a typi- 
cal equal-mass, nonspinning MBHB source of this mass 
and redshift within a single LSST field several months in 
advance of merger. This will allow multiple exposures of 
the target area to be made throughout the event in an 
effort to detect an optical transient. 

If a transient signal is not bright enough to be visible 
with a survey telescope such as LSST, then a more sen- 
sitive instrument with a narrower field of view may be 
required. In Fig. 1 we plot the 2.2 arcmin x 4.4 arcmin 
FOV of the NIRCam instrument on the James Webb 
Space Telescope (JWST) [36] as an example of such an 
instrument. Our most complete estimate of LISA’s lo- 
calization capabilities (IMR waveform, orbiting detector 
with Synthetic LISA response, all TDI channels) has 
a median latitude error of 2.6 arcmin. This suggests that 
a single JWST exposure might be sufficient to cover the 
LISA error box. In reality, a few exposures may be re- 
quired to account for the shape of the LISA error box, 
its orientation with respect to the JWST FOV, and any 
operational constraints of JWST. The results in Fig. 1 
suggest that unlike the situation with LSST, in which a 
number of pre-merger images could be obtained, a co- 
ordinated observation with JWST (or any other instru- 
ment with a similar FOV) will have to be a follow-up, 
as LISA requires the complete IMR waveform to provide 
precise localization information. For example, at ISCO, 
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FIG. 2: Latitude uncertainty, a 13 , for a 2 x 10 e M© system at 2 = 1. The center line (with individual time samples marked) 
denotes the median error, while the shaded region depicts the middle three quintiles (20 th — 80 fh percentile) of errors. Both 
plots correspond to IMR signals with an orbiting LISA constellation. However, Panel (a) shows results for just the A and E 
TDI channels, while Panel (b) also includes the T channel. 




FIG. 3: Same as Fig. 2, but for uncertainty in longitude, a x . 


the best estimate for median latitude error exceeds the 
JWST FOV by roughly half an order of magnitude. 

All of the curves in Fig. 1 represent the median LISA 
latitude error for an ensemble of MBHB mergers. Fig- 
ures 2 and 3 also show the distribution of both lati- 
tude and longitude errors, based on estimates with the 
complete IMR. waveform and LISA’s orbital modulation. 
Separate panels show results including just A and E and 
results including all three TDI channels. In each panel, 
the solid lines represent the median values (equivalent to 
what is plotted in Fig. 1), while the shaded area shows 
the region spanned by the middle three quintiles of the 
distribution. The left panels {A, E only) show that the 
final LISA error box will lie within the FOV of JWST for 
systems in the 80 ,ft percentile of sky localization preci- 
sion. That is, both latitude and longitude are measured 
to within ~ 4.4 arcmin or better for the top 20% of sys- 
tems. When T is included, the situation is better still, 
with a full 50% of binaries in our ensemble fitting within 
the JWST FOV 1 . These figures show just how important 
T is for computing final localization errors. Not only do 


1 These assertions are somewhat crude: The plot does not actually 
show whether individual binaries meet the criterion. It is possi- 
ble that some systems are well localized in latitude but poorly 
localized in longitude, and vice-versa. A more careful check of 
the raw data verifies the numbers in the text. 


the medians drop for both parameters, but the width of 
the distribution decreases as well. Although a 5-link and 
6-link LISA are indistinguishable at early times, com- 
plete functionality can have important effects for the full 
signal. 


For completeness, in Fig. 4 we show errors in lumi- 
nosity distance, a Dc -/Di, for IMR signals with orbital 
modulation. Distance is the third axis of the 3D local- 
ization “voxel.” Given an assumed cosmological model, a 
measurement of distance is equivalent to a measurement 
of redshift. This approximate redshift can help decide be- 
tween potential host galaxies located within the 2D error 
box. We see that the median distance error drops below a 
percent within the last day, reaching a minimum less than 
0.1%. Unfortunately, this “intrinsic LISA error” is mis- 
leading. Weak leasing from matter between the source 
and detector will slightly magnify or demagnify the GW 
signal in an unknown way, adding significant error to the 
distance measurement. The line marked “WL” on the 
figure is an estimate [341 of <j Dl /Di ~ 4.4% from weak 
lensing at 2 = 1. The intrinsic LISA error crosses the 
WL error over a month before merger; shortly after that 
time, the information on provided by LISA becomes 
irrelevant. 
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FIG. 4: Same as Fig. 2, but for fractional luminosity distance uncertainty, a° L /D l. An estimate of the weak tensing limit 
( WL) [34] at z = 1 is shown for reference. Figs. 2-4 collectively show the evolution of the 3D error voxel. 


IV. CONTRIBUTION OF MERGERS TO 
PARAMETER ESTIMATION 

The results presented in the preceding section clearly 
demonstrate that observing the MR portion of a MBHB 
signal can provide a significant amount of sky localization 
information, despite the fact that there is little or no 
modulation imposed by the LISA orbits. In this section, 
we explore the nature of this information. 

One question of interest is whether the information 
provided by the MR signal is tightly coupled to just one 
or two of the elements in A“ or if it generally improves 
the knowledge of all parameters. Figure 5 shows the time 
evolution of the median eigenvalues of the Fisher ma- 
trix for the ensemble of systems studied in Section III. 
(Throughout this section, we always use the complete 
IMR waveform with orbital modulation and all three TDI 
channels.) At each time step, T a i, was calculated using 
the A, E, and T observables for every case in the ensem- 
ble using the procedure described in Section II. These 
matrices were used to compute the eigenvalues for every 
case in the ensemble, and the median values are plotted. 
A similar procedure was used to compute the median 
eigenvectors corresponding to these eigenvalues. 

As (8) indicates, the Fisher matrix is a measure of 
inverse error or “information.” The increasing trends in 
the curves in Fig. 5 represent information accumulating 
in the LISA observables as the signal evolves, just as 
the decreasing trends in Figs. 1-4 represent a decrease 
in measurement error. The largest eigenvalue in Fig. 5 
(solid curve) corresponds to an eigenvector that mostly 
consists of one parameter, In M. The majority of the 
information in this eigenvalue accumulates early in the 
inspiral phase of the signal. From ISCO through the end 
of the signal, the eigenvalue roughly doubles, with most of 
the increase occurring very near the time of coalescence. 
This is consistent with previous results that show that 
hi M is primarily determined during the inspiral phase 
137] and that the MR signal provides little additional 
information [15]. By comparison, the increase in SNR 2 , 
which one might naievly expect to be proportional to 
Fisher matrix eigenvalues, is ~ 100 over this period. 

The band of eigenvalues in the middle of Fig. 5 (dasiied 



FIG. 5: Evolution of the median Fisher matrix eigenvalues 
with time for the inspiral, merger, and ringdown of an ensem- 
ble of equal-mass, nonspinning MBHBs with M a = 2x 10®Mg 
at 3=1 observed by an orbiting, six-link LISA via the A, E, 
and T TDI channels. 


curves) represent a mix of the five angular parameters in 
A° and In Di. The increase in these eigenvalues from 
ISCO through ringdown is much more significant than 
that for In M, ranging from two to three orders of mag- 
nitude. 

The smallest eigenvalue (dotted curve) corresponds to 
an eigenvector that mostly consists of t c . The small size 
of this eigenvalue is partly due to the intrinsic scale of the 
parameter, which was the only one in A° to be expressed 
in physical units (seconds). It increases by a factor of 
~ 30 between ISCO and the end of ringdown. 

The fact that all of the eigenvalues increase post-ISCO 
indicates that the MR portion of the waveform, coupled 
with the high-frequency response of the LISA instrument, 
provides information in all directions of parameter space 
(although the information content in the InM direction 
is a bit more modest than in the other directions). This 
behavior contrasts with some previous speculations that 
the primary effect of the MR signal would be to constrain 
t c via the sharp peak that occurs in the waveform near co- 
alescence. In that scenario, any improvement in other pa- 
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rameters would arise via reductions of correlations with 
t c . Figure 6 investigates this specific mechanism in more 
detail. For each parameter plotted, there are three sets of 
histograms representing parameter errors for each case in 
the ensemble. The first (solid) includes the entire signal. 
The second (dashed) is obtained from signals truncated 
in time at ISCO. The third set (dotted) is an estimate 
of how the ISCO histograms would be affected should an 
exact measurement of t c be provided. This was obtained 
by removing the row and column corresponding to t c in 
the Fisher matrix prior to applying (8). 

With the exception of latitude 3, the dotted his- 
tograms track tlie dashed histograms and not the solid 
ones. This indicates that the general improvement in 
parameters that occurs when the entire MBHB signal is 
considered is not due to the MR, signal providing an ac- 
curate determination of coalescence time. 

In the case of latitude, the dotted histogram lies 
halfway (in log space) between the dashed and solid his- 
tograms. A correlation between t c and sky position might 
be expected because of the definition of t c used in our 
study. We measure the time of coalescence at the solar 
system barycenter, while it is more naturally measured 
as the time the peak hits the detector. The difference in 
arrival times depends on the sky position. However, we 
would expect this correlation to affect both ecliptic lati- 
tude and ecliptic longitude, so it is not clear why latitude 
responds differently. 

V. SUMMARY AND CONCLUSIONS 

We have presented the first estimates of LISA’s ability 
to measure the source parameters of nonspinning, equal- 
mass MBHBs that include a complete inspiral-merger- 
ringdown waveform, realistic LISA orbits, and the full 
high-frequency response of the LISA instrument. For sys- 
tems with M 0 = 2 x 10 6 Mq at z — 1, we find that LISA 
can locate the source on the sky with an error of ~ 3 
arcminutes in the median case at the end of the MBHB 
event. We have also studied how the sky localization in- 
formation for these systems evolves with time and find 
that, at ISCO. the sky localization errors are typically 
tens of arcminutes. 

Inspired by the potential for multi-messenger astron- 
omy enabled by combined observations of MBHB mergers 
with gravitational-wave detectors and electromagnetic 
telescopes, we have compared the source-localization ca- 
pabilities of LISA with the FOVs of two future telescopes: 
LSST and JWST. We find that most systems are located 
within the LSST FOV several months prior to coales- 
cence. Localization within the smaller FOV of JWST 
requires information from the complete waveform. 

While the results in this paper focus on a particular 
instrument (LISA) and a particular source class (equal- 
mass, nonspinning, M a = 2 x IQ 6 M®, z = 1) they 
can give some insight to other scenarios. For example, 


proposals for lower-cost variants of LISA often involve 
a constellation with shorter arm lengths [38]. This has 
the general effect of shifting the response of the instru- 
ment to higher Fourier frequencies. We have asserted 
that the improvement in parameter estimation that we 
observe post-ISCO is due in part to the complexity of 
the high-frequency response of the instrument. If this 
is true, shortening the arms could have a significant im- 
pact on the ability to localize systems in the mass range 
studied. Specifically, the improvement in sky localization 
from observing the merger-ringdown of a system with 
M a ~ 10® M© may not be as dramatic as that presented 
in Section III. At the same time, the decrease in strain 
sensitivity associated with shorter arms would reduce the 
duration of the inspiral signal that is observed, further 
degrading measurement performance for these systems. 
For systems with larger M 0 , the situation would be worse. 

For lower-mass systems, shortening the arms (if the in- 
terferometric measurement system is not simultaneously 
degraded) may actually increase the benefits of observ- 
ing the merger by increasing the SNR for the merger- 
ringdown portion of the signal. Since these lighter sys- 
tems have higher merger frequencies, they will still ex- 
perience a complex response in the detector, even with 
shorter arms. Determining whether this effect can over- 
come the lower intrinsic SNR of less massive systems will 
require a more detailed calculation. 

Finally, we note that while the addition of merger + 
ringdown and a complete detector response is an impor- 
tant step towards realizing the full potential of LISA and 
other space-based gravitational-wave detectors, we still 
do not have a complete picture. Other potential sources 
of information in the MBHB waveform, particularly spin 
and tiie associated precession effects, have been demon- 
strated to significantly improve LISA’s localization abil- 
ity in the inspiral-only case [7]. It remains to be seen 
what impact combining spin precession and merger will 
have on LISA sky localization estimates. 
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FIG. 6: Comparison of sky localization accuracy including the complete signal (solid), only the inspiral (dashed), and the 
inspiral combined with the assumption of perfect knowledge of t c (dotted). All histogram bins are normalized by the total 
number of cases and are expressed as percentages. 
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